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dimension  of  U-shaped  configuration  defined  in  Figure  8 
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average  value  of  D'  during  application  of  load  increment 
value  of  D'  at  end  of  i-th  element 
modulus  of  elasticity 

effective  elastic  modulus  defined  by  Equation  (42) 
functional  relation  of  Equation  (45) 
shear  modulus 

torsional  rigidity  per  unit  length  of  cross  section  of  axial  strip 
in  convoluted  cylindrical  structure 

thickness 

index 
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defined  by  Equation  (80) 
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length  of  cantilever  beom  of  Figure  1;  width  of  oxial  strip  defined 
in  Figure  25 
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moment  per  unit  length 
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differential  pressure  across  convoluted  cylindrical  structure 

initiol  radius  of  convolute 

rodius  of  cone 

average  radius  of  convoluted  cylindrical  structure 
coordinate  of  length  along  convoluted  curve 
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NOMENCLATURE  (cont*d) 


load  por  unit  length 

2 

reference  load  per  unit  length,  D/A* 
overoge  value  of  T  during  application  of  load  increment 
value  of  T  after  application  of  n-th  load  increment 
n-*th  load  increment 

elastic  stroin  energy  of  bending  per  unit  width  of  cantilever  beom 
of  Figure  1 

volue  of  U  after  applicotion  of  n*th  load  increment 

bending  energy  in  convoluted  cylindrical  structure 

bending  energy  required  to  form  convoluted  cylindricol  structure 
from  straight  convoluted  strip 

twisting  energy  in  convoluted  cylindrical  structure 

total  energy  in  convoluted  cylindrical  structure 

axial  coordinate  of  convoluted  cylindrical  structure 

Cartesian  coordinotes  with  x*axis  tongent  to  convoluted  curve  for 
given  load;  coordinates  of  Figure  4  used  in  plastic  bendir>g  orKilysis 

values  of  x,y  corresponding  to  initial  undeformed  structure 

values  of  x,y  ofter  application  of  load  increment  (Figure  2) 

moment  orm  in  two-dimensiorKil  structure 

value  of  z  at  end  of  i-th  element 

average  volue  of  Z|  during  applicotion  of  lood  increment 

maximum  value  of  z  for  a  given  lood 

value  of  ofter  opplicotion  of  n«*th  lood  increment 
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cone  half-angle 

parameter  defined  by  Equations  (19)  or  (26) 

porameter  defined  by  Equation  (27) 

constants  defined  by  Equations  (70) 

end  deflection  of  elemental  cantilever  beam 

end  slope  of  elemental  cantilever  beam 

deflection  of  end  of  i-th  element  relative  to  (i  -  1)th  element 

slope  of  end  of  i-th  element  relative  to  (i  -  1)th  element 

plastic  bending  strain 

"plastic  stroin"  at  proportional  limit  defined  by  Equations  (48)  and  (49) 

slope  of  tangent  to  two-dimensional  structure 

initial  value  of  0  corresponding  to  unloaded  structure 

value  of  0  at  end  of  i-th  element 

average  value  of  0.  during  opplication  of  load  increment 
chonge  in  slope,  0  *  9^ 
porometer  defined  by  Equation  (41) 

Poisson's  ratio 

horizontal  and  vertical  deflections,  respectively,  of  end  of 
cantilever  beom  of  Figure  1 

horizontal  and  verticol  deflections,  respectively,  of  end  of  i-th 
element  of  Figure  1c,  relotive  to  (i  -  1)th  element 

incremental  change  in  ^  during  application  of  n-th  load  increment 


NOMENCLATURE  (cont'd) 


bending  stress 
maximum  bending  stress 
proportiorKil  limit  stress 
yield  strength 

angle  defined  in  Figure  25;  change  in  slope  of  S*shaped  strip  of 
convoluted  cylindrical  structure 

maximum  value  of  for  S*shaped  strip  of  convoluted  cylindrical 
structure 

minimum  and  maximum  values  of  <P  ,  respectively,  defined  in 
Figure  25 

material  efficiency  factor  defined  by  Equation  (56) 
value  of  after  application  of  n*th  load  increment 
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FINAL  REPORT  ON  "RESEARCH  ON  ENERGY  ABSORBING  STRUCTURES 


Theoretical  and  experimental  studies  of  convoluted  type  structures  ore 
described.  An  incremental  method  for  predicting  the  nonlineor  lorge-deformotion 
behavior  of  two*dimensional  convoluted  strips  is  shown.  The  method  is  opplicable 
for  elostic  as  well  as  plastic  stress-stroin  behavior  of  the  material.  The  results  of 
the  two'^imensionol  artolysis  are  extended  to  an  analysis  of  o  circumferentially* 
convoluted  cylinder  which  deforms  into  o  conical  shape.  Tests  on  typical  two*  and 
three-dimensionol  convoluted  metol  specimens  are  described  and  correlations  of 
theoretical  predictions  with  experimentol  results  ore  shown. 

II.  INTRODUCTION 

Flexible  metal  structures  which  can  grouly  change  their  shape  by  the 
obeorption  or  release  of  energy  to  their  environment  offer  an  attractive  potentiol 
for  aerospoce  applicotions.  In  particular,  it  would  be  desirable  to  be  able  to  utilize 
o  passive  metol  structure  of  preset  convoluted  geometry  which,  when  acted  upon  by 
extemol  environmental  forces,  would  vary  its  geometry  in  o  predetermined  manner. 
These  forces  can  orise  from  numerous  sources;  e.g.,  oerodynamic  pressures,  thermal 
stresses,  voriations  in  materiol  properties  due  to  temperature  chortges,  etc.  In  view 
of  the  limited  knowledge  avoiloble  on  this  type  of  structure  o  theoreticol  end  experi¬ 
mental  program  is  being  carried  out  to  develop  methods  suitoble  for  their  onolysls. 


III.  SCOPE  OF  WORK 


The  scope  of  the  work  performed  under  this  study  is  summarized  ot  follows: 

(1)  An  incremental  technique  wos  developed  for  predicting  the  large- 
deflection  lood-deformation  behavior  of  two-dimensioiKil  convoluted  metal  strips 
of  arbitrary  shape  which  exhibit  elostic  and  plostic  material  behavior. 

(2)  The  load-deformation  behavior  of  two-dimensional  convoluted 
strips  loaded  within  the  elastic  range  wos  computed  for  several  typicol  configura¬ 
tions  ond  is  presented  in  nondimensional  form. 

(3)  The  load-deformation  behavior  of  two-dimensionol  convoluted 
strips  of  semi-circular-orc  configurotion  was  conqsuted  for  several  conditions  of 
yielding  and  plastic  stress-strain  behavior  ond  is  presented  in  nondimensionol  form. 

(4)  The  results  of  the  two-dimensional  analyses  were  exter>ded  to  on 
analysis  of  a  three-uimensional  circumferentially-convoluted  cylindrical  structure 
which  deforms  into  a  conicol  shape. 

(5)  Tests  were  performed  on  typical  two-  ortd  three-dimensionol 
structures  and  the  results  were  correlated  with  the  theoretical  predictions. 

IV.  TWO-DIMENSIONAL  ANALYSIS 

The  large-deflection  onolysis  of  convoluted  metal  structures  is  severely 
complicoted  by  the  fact  that  the  lorge  deflections  give  rase  to  nonlineorities  in  the 
load-deflection  behavior.  Even  when  the  moteriol  stress-stroin  behovior  is  elostic, 
geometrical  nonliraearities  ore  present.  If,  in  oddition,  the  moteriol  is  deformed 
into  the  plostic  ronge,  nonlineorities  result  in  the  streu-stroin  behavior  ond  the 


probUm  It  further  compMcaftd.  Because  of  these  complications,  onolytical 
closed-form  solutions  con  only  be  obtained  for  the  simplest  cases,  and  severe 
assumptions  and  approximations  are  usually  required. 

In  order  to  arrive  at  more  general  and  useful  results,  an  incremental 
approach  hot  been  taken  in  the  analysis  of  two-dimensional  convoluted  strips. 
Solutions  ore  developed  in  increments  of  lood,  and  the  incremental  deformations 
resulting  In  each  lood  increment  are  determined  by  dividing  the  structure  into 
incrementol  segments.  This  technique  is  described  in  detail  in  the  following 
sections. 

A.  Development  of  Incrementol  Method 

Consider  the  simple  two-dimensional  strip  shown  in  Figure  la, 
which  is  deformed  by  loads  applied  to  its  ends.  As  the  load  increases,  deformation 
results  from  changes  in  the  bending  moment  distribution  and,  if  the  material  is  loaded 
into  the  plastic  range,  changes  in  the  effective  flexurol  rigidity  of  the  structure. 
Since  large  deformations  are  being  considered,  changes  in  the  bending  moment  ot 
any  point  are  due  to  changes  in  the  load  as  well  os  changes  in  the  effective-moment 
orm.  Because  of  symmetry  it  is  necessary  to  analyze  only  one  typical  section  of  the 
strip,  such  as  the  curved  contilevered  beam  of  Figure  1b.  At  a  point  "a”  the  change 
in  slope  is  zero  arul  at  point  ”b"  the  moment  is  zero  since  the  moment  arm  vanishes. 
As  the  structure  deforms,  the  edge  load  T  (load/unit  length)  remoins  horizontal  so 
that  the  berKling  moment  ot  any  point  is  the  product  of  T  ond  the  moment  orm,  z. 

The  basic  problem  is  to  describe  the  deformation  of  a  convoluted 
strip  of  arbitrary  shape  os  the  lood  increoses  from  zero  to  some  arbitrary  volue.  For 
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FI6.  I  GEOMETRIC  PA^RAMETERS  AND  COORDiNIkTVE 
USED  IN  THE  A.NA»UVSI»  OF 
TWO  ■  Ol  MEN&ION«kU  STRIPS 


(d) 


this  purpose  tht  following  ossomptlons  oro  modo: 

1 .  Tho  thlcknoss  of  tho  motoriol  is  small  onough  in  rolotlon  to  tho 
othor  dimonsions  that  clouicol  bonding  thoory  con  bo  usod  to  roloto  chongo  in 
curvaturo  to  bonding  momont. 

2.  Tho  offoch  of  normal  forcos  acting  in  tho  tangont  piano  at  ony 
point  con  bo  ignored. 

3.  Motor iol  stross-stroin  behavior  in  both  tho  olostic  ond  plastic  range 

occurs. 

4.  Tho  structure  is  doformod  in  o  stoto  of  piano  strain  bonding  and 
ond* effects  con  bo  ignored. 

With  tho  above  assumptions  tho  deformation  bohovior  is  determined  by  o 
stop*by*stop  procedure  in  which  on  incremental  lood  is  oppliod  and  tho  incremental 
deformation  is  ostimotod  from  tho  inital  shape  ond  the  oppliod  load.  Since  tho 
moments  and  (in  tho  plastic  range)  tho  effective  floxurol  rigidity  chongo  during 
application  of  tho  lood,  overage  values  ore  usod  which  con  bo  determined  by  triol 
and  error.  In  order  to  determine  the  incrementol  deformotions  a  second  incrementol 
procedure  is  utilized  in  which  the  curved  contilever  boom  of  Figure  1b  is  opproximoted 
by  straight  segments,  os  shown  in  Figure  1c.  The  onolyticol  development  of  both 
incrementol  procedures  is  given  below. 

The  bending  equotion  con  be  expressed  by 
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wh«r«  0  is  th#  slop«,  0^  is  the  slope  for  the  initial  undeformed  shope,  s  is  the 
coordirKite  of  length  along  the  curved  section,  and  D'  is  an  effective  flexural 
rigidity. 

The  moment  M  is  given  by 


(2) 


For  elastic  material  behavior. 


O  ^  — - :  , 


(3) 


where  E  is  Young's  modulus,  h  is  the  material  thickness,  and  V  is  Poisson's 
ratio. 

Equation  (1)  expresses  the  change  in  curvature  due  to  the  moment  M. 

The  curvature  dO/ ds  can  also  be  expressed  in  terms  of  Cartesian  coordinates  x,y 
by  the  relation* 

c/^y 

Since  the  location  of  the  x,y-axes  ii  arbitrary,  it  is  useful  to  consider  a  "moving" 
set  of  coordinate  axes  with  the  origin  on  the  curved  beam  and  the  x-axis  tangent 


*  it  is  conventent  to  select  the  coordinate  axes  so  that  a  positive  change  of  curvature 
in  the  x,y  coordinates  corresportds  to  a  decrease  in  the  original  curvature  of  the 
convoluted  shape. 
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to  tho  boom  ot  tho  point  $•  Tho  first  dorivotivo  form  in  Equotion  (4)  is  thon  zoro 
and  Equations  (1)  and  (2)  moy  bo  writton 


^  Us 

<i/x*  ^  ^ 


whoro  tho  subscript  o  donotos  tho  initial  undoformod  shape.  Equation  (5)  expresses 
tho  curvature  rolotion  at  a  particular  point  s  since  the  coordinate  axes  of  x,y  are 
now  functions  of  s  and  T,  and  those  of  x^,y^  correspond  to  T  =  0.  However, 
if  Equotion  (5)  is  extended  to  a  smoll  finite  range  of  x  about  x  ~  0  for  fixed 
coordinate  axes  of  x,y  and  x^,y^,  the  error  will  be  small  since  the  first  derivative 
term  in  Equation  (4)  will  be  small  compared  with  unity.  This  is  the  basis  for  the 
present  incremental  procedure. 

Consider,  first,  an  incrementol  change  in  the  load  and  express  the 
resulting  change  in  curvature  from  Equation  (5).  This  result  can  be  written 


D'  MT 


(6) 


providad  the  load  increment  is  sufficiently  small.  The  term  on  the  left  of 
Equation  (6)  can  be  interpreted,  with  the  aid  of  Figure  2,  os  follows:  Let  x,y 
represent  the  coordinates  at  point  s  which  correspond  to  the  load  T,  and  let 
x',y'  correspond  to  T  +  Zi  T,  os  shown  in  Figure  2a.  By  superposing  the 
coordinate  axes,  os  itKiicated  in  Figure  2b,  the  incremental  change  in  curvature 
is  seen  to  be 


ri6. 


(a) 


(  b) 


dtOMKTItIC  INTBItPRKT^^TION  OF  LCFT  SIDE 
or  EQUiKTION  C6) 


(7) 


•'»*  ■  <'«’  ■^«*  "izr^  =  't;^  •> 

wh«r»  is  the  incremental  deflection  normal  to  the  tangent  at  the  point  $. 
Equations  (6)  and  (7)  can  be  written 

*  O'  *  (8) 


where 


AM  s  r*  t  T 


Jr 


-  Ti 

t>* 


olT 


)^T 


(9) 


The  quantity  is  an  effective  incrementol  moment  which  also  takes  into 

occount  variotions  in  D*  which  result  if  the  material  is  in  the  plastic  range. 

Equations  (8)  and  (9)  form  the  basis  for  the  incrementol  technique  by  which 
the  incremental  deformation  is  determined  for  each  lood  increment.  As  mentioned 
eorlier,  the  cantilever  beam  of  Figure  1b  is  opproximated  by  a  series  of  straight 
elements,  shown  in  Figure  1c.  If  the  increments  of  length  and  load  are  smoll 
enough,  each  element  con  be  treated  os  o  contilever  beam  with  linearly  varying 
moment  ond,  in  the  plastic  range,  linearly  varying  flexurol  rigidity.  The  incremental 
changes  in  lateral  deflection  and  slope  of  the  end  of  eoch  element  relotive  to  its 
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boM  0T9  dateri^inad  from  Equations  (8)  and  (9),  The  gross  deflections  are  then 
determined  by  accumulating  rotations  and  deflections  of  each  element  starting  from 
the  base  of  the  original  contilever  beom. 


Let  the  relative  incremental  slope  and  bending  deflection  for  the  i-th 
element  of  Figure  1c  be  denoted  by  and  respectively*  Positive  values 
for  and  correspond  to  a  clockwise  rotation  of  the  element,  as  shown  in 
Figure  3.  The  total  lateral  deflection  of  the  end  of  the  i-th  element  relative  to 
the  end  of  the  (I  -  1  )th  element  results  from  and  the  rigid  body  rotation 

produced  by  the  bending  of  all  the  preceding  elements,  and  is  given  by 


totol  lateral  deflection  of  end  of  i*th  element 
relative  to  end  of  (i  -  1)th  element 


-  2  a 


where  is  the  length  of  the  i*th  element.  The  average  slope  of  the  i-th  element 
6*  during  application  of  the  load  increment  can  be  approximated  by 


wh*ra  0|  is  the  slope  before  application  of  the  load  increment. 


The  horizontal  and  vertical  incremental  defections  of  the  end  of  the  i-th 
element  relative  to  the  (i  -  1)th  element  ore,  from  Equations  (10)  and  (11), 


*  In  terms  of  tlie  previous  notation,  m  ^  9  4^1  . 


-9- 


With  the  fixed  end  boundary  conditions  for  the  original  cantilever  beam  denoted  by 


the  incremental  deflections  of  the  end  of  the  i-th  element  con  be  written 

A?,  =  2 

fi-l 

Equations  (12)  -  (14)  constitute  the  relations  from  which  the  gross  deflections  of 
the  original  cantilever  beam  can  be  computed  for  eoch  load  increment,  provided 
the  quantities  and  are  known.  The  determination  of  £•  and  for 
elastic  and  plastic  material  behavior  is  described  in  the  following  sections. 


B.  Incrementol  Slopes  and  Deflections  for  Elostic  Moterial  Behavior 
For  elastic  behavior  Equation  (3)  applies  so  that  dD'/ dT  is  zero 


and  Equations  (8)  and  (9)  become 


(15) 


where 

m  -  (f+r 


(16) 


The  quontitles  ^  and  j !  can  be  determir>ed  froAn  these  two  equotions.  Since 


finitt  variations  occur  in  z,  T,  and  dz/dT  during  application  of  the  load  increment, 
averoge  values  of  these  quantities  are  used  to  compute  AM. 

As  was  mentioned  earlier,  if  the  increments  of  length  are  small  enough  the 
voriation  in  AM  over  each  elemental  beam  can  be  assumed  linear.  Thus,  let 

where 


(Zr 


(18) 


and  the  bars  denote  average  quantities  during  application  of  the  load  increment.  The 
subscript  i  on  M  and  z  denotes  a  quantity  evaluated  at  the  end  of  the  i-th 
element,  and  oi.  is  defined  by 


- 


AM,.,  -  4/^.- 


(19) 


Substitution  of  Equation  (17)  into  Equation  (15),  and  integration  of  the  result  with 
the  boundary  conditions. 


Ay(o)  ~ 


(20) 


yields 


- 


c/k 


«■)  ■=  ('/-*) 


(21) 
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and 


tyUi) 


m--,  <■' 

Zb 


(1- 


(22) 


Equations  (21)  and  (22)  with  Equations  (18)  and  (19)  constitute  the  required 
relations. 

For  computational  purposes  it  is  convenient  to  introduce  the  reference 
load  and  moment. 


(23) 


and  to  express  the  above  equations  in  nondimensional  form.  Here  A,  is  a  typical 
dimension  of  the  convolute,  such  os  on  initial  radius  of  curvature.  Equations  (21) 
and  (22)  become 


/ 


(24) 


ond 


k 

4. 


(25) 
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C.  Incrtmantal  Slopes  and  Deflections  for  Plastic  Material  Behovior 

In  rhe  pigstic  case  both  AM'  and  D'  vary  over  the  length  of  the 
elenoent.  As  before,  if  the  element  length  is  sufficiently  small,  both  of  these 
quontities  can  be  assumed  to  vary  linearly.  Thus,  let 


-4/V/ 


and 


J 


(26) 


(27) 


so  that 


and 


(  /- 


(28) 


(29) 


Substitution  of  Equations  (28)  and  (29)  into  Equation  (8),  and  integnation  of  the 
result  with  the  boundary  conditions  of  Equations  (20),  gives 


(30) 


(31) 


Ai  before,  it  is  convenient  to  use  Equation  (23)  and  to  write  Equations  (30)  and 
(31)  in  the  nondimensionol  form. 


^  /o 


(32) 


A 


(** 


of-^ 


(33) 


If  the  volue  of  is  small  these  equations  are  difficult  to  evaluate  due 
to  small  differences  between  large  numbers.  For  this  case  the  following  series 
expansions  are  useful: 


(34) 


—  -  ^ )  r  /  »  (?•  s’- 


(35) 


D.  Plostic  Bending  Re  lotions 

In  the  plostic  onalysis  described  above  a  knowledge  of  the  moment- 
curvature  relation  is  required  in  order  to  calculate  the  effective  flexural  rigidity 
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and  the  porameter  Becouse  of  the  numerical  procedure  involved  this  relation 

con  be  determined  from  the  actual  $tre$$*strain  behavior  of  the  material  and  utilized 
in  graphical  or  tabular  form.  However,  in  order  to  illustrate  the  method,  on  ideal 
elastic-plastic  moterial  behavior  is  first  assumed,  for  which  an  analytical  solution 
can  be  obtained  in  closed  form.  The  following  analysis  is  based  on  that  described 
in  Reference  A. 


1 .  Ideal  Elostic-Plastic  Behavior 

Consider  the  bending  of  a  uniform  sheet  of  nonhardening 
material  under  conditions  of  plane  strain.  The  material  is  assumed  to  have  the 
stress-strain  behavior  of  Figure  4a  in  both  tension  and  compression  and  is  assumed 
to  yield  in  accordance  with  the  Tresca  law.  It  is  further  assumed  that  the  radius 
of  curvature  is  so  lorge  in  relation  to  the  thickness  that  the  induced  transverse 
stresses  in  the  direction  normal  to  the  sheet  can  be  neglected. 

The  coordinate  axes  are  chosen  as  shown  in  Figure  4b  with 
the  z-oxis  in  the  direction  in  which  strain  is  prevented.  When  the  moment  per  unit 
width  is  greater  than  the  yield  value, 

^  6 


elastic  and  plastic  zones  are  produced,  as  shown  in  Figure  4b.  The  moment  is 


expressed  by 

2E 


o 


c 


zee 


3 


SO-V‘3 


ds 


(37) 
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Figur*  4.  FaraiMtan  ond  CoordinatM  Ihcd  in  Idaolind  PIcnHc  Banding  Analysli 


I 


-17- 


ond  Hi«  parameter  c  is  given  by 


( l~  v'’)  op 

Solution  of  d(0  *  Q^)/ ds  from  Equations  (37)  and  (38)  gives  the  result, 


dlC6-6,^ 


(38) 


(39) 


Thu  result  with  Equations  (1),  (3),  and  (36)  gives  for  the  effective  flexural 
rigidity. 


My  ^  . 


(40) 


With  this  material  behavior  one  additional  structurol  parameter  is  required 
in  order  to  describe  the  load  or  deflection  value  at  which  the  material  first  yields. 
For  this  purpose  it  is  convenient  to  introduce  the  parameter  ^  ,  defined  by 


(41) 


where  is  the  proportional  limit,  equal  to  ^  for  this  case,  and  Ey  is 

defined  by 


(42) 
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Equation  (40)  with  Equotions  (3)  and  (23)  can  then  be  written. 


Equations  (40),  or  (41)  and  (43),  constitute  the  required  relations.  A  treatment  of 
the  problem  for  more  general  stress-strain  behavior  is  given  in  the  following  section. 

2.  Arbitrary  Stress-Stroin  Behovior 

In  the  previous  section  relations  were  derived  for  plane 
strain  bending  of  a  wide  sheet  of  ideal  nonhardening  elastic-plastic  material.  The 
present  ar>alysis,  based  on  that  of  Reference  B,  extends  the  treatment  to  arbitrary 
stress-strain  behavior.  As  in  the  former  analysis,  the  following  assumptions  are  made: 

(i)  The  material  has  the  some  stress-strain  behavior  in 
tension  and  compression* 

(ii)  Transverse  stresses  in  the  direction  normal  to  the 
sheet  con  be  neglected. 

(iii)  The  neutral  axis  remains  coincident  with  the  midplane 
of  the  sheet. 

(iv)  The  total  strain  compor^nt  in  the  restrained  direction 
is  zero. 

In  addition  to  the  above,  the  usual  lows  of  plasticity  are 
imposed,  which  include  the  constancy  of  volume  condition,  the  von  Mises  condition, 
and  a  criterion  of  plastic  flow.  For  the  letter  the.  "moximum  shear  stress  vs. 
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numtrlcally  largest  principal  strain"  critarion  is  used.*  With  these  conditions  the 
result  of  Reference  B  may  be  written  as 


where  CT  is  the  bending  stress  and  Cp  is  the  plastic  bending  strain. 

In  addition  to  Equation  (44)  there  exists  a  stress-strain 

relation  between  and  0^  , 

P 


-  j 


(45) 


which  can  be  determined  from  a  tensile  test.  The  present  flow  criterion  is  such  that 
the  functional  relotion  of  Equation  (45)  is  the  same  as  that  between  oxial  stress  and 
plastic  strain  in  a  tensile  test. 

Equations  (44)  and  (45),  with  the  expression  for  moment 


per  unit  width. 


2  I  (Tycfy 
o 


(46) 


can  b«  uMd  to  compute  the  moment-curvature  relation.  For  a  particular  value  of 
d60/dt,  values  of  CT  can  be  selected  for  which  the  corresponding  values  of  y 
can  be  determined  from  Equations  (44)  and  (45).  Equation  (46)  can  then  be  used 
to  compute  the  corresponding  moment. 


*  Although  the  von  Mises-Hencky  criterion  is  generally  considered  to  be  slightly 
better,  the  added  complexity  Is  not  worranted  In  this  oiwlysis. 
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If,  instead  of  Auumption  (iv),  which  is  the  usual  condition 
for  plane  strain,  it  is  auumed  that  both  the  elastic  and  plastic  components  of  strain 
ore  zero,  the  result  is 

~  J  (47) 

which  it  considerably  simpler  than  Equation  (44).  A  comparison  of  Equations  (44) 
and  (47)  is  given  below. 

For  many  structural  metals  the  stress-strain  relation  of 
Equation  (45)  can  be  expressed  in  the  form 

(T "  .  (48) 

With  the  addition  of  the  elastic  strain,  OVE,  the  expression  for  total  strain  becomes 
identical  to  the  well-known  Ramberg-Osgood  relation  (Reference  C)  by  proper 
definition  of  the  constant  A.  It  is  convenient  to  express  Equation  (48)  in  terms  of 
the  proportional  limit  stress  and  some  corresponding  plastic  strain  This 

result  can  be  written  as 


Although,  strictly  speaking,  there  is  no  plastic  strain  at  the  proportional  limit, 
Lubahn  and  Felgar*  have  suggested  20  x  10  ^  as  a  suitable  value  since  this  is 
close  to  the  minimum  value  that  can  be  detected  experimentally. 


*  Reference  B,  p.  V1 1 
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In  Reference  B  a  comparison  of  Equations  (44)  and  (47)  is 
given  for  a  porticular  case  of  bending,  using  the  expression  of  Equotion  (48)  with 
n  =  5*  Further  comparisons  are  shown  in  Figure  5  for  n  =  3,  5,  and  10.  For  each 
of  these  cases  the  maximum  plastic  bending  strain  is  half  of  the  total  strain.  From 
these  results  it  appears  that  Equation  (47)  provides  a  good  approximation  to 
Equation  (44). 

With  the  stress-strain  relation  of  Equation  (49)  analytical 
expressions  con  be  determined  from  Equations  (46)  and  (47)  which  relate  di&O/ds 
and  M.  Substitution  of  Equation  (49)  in  Equation  (47)  gives 


(50) 


2 

where  E^  =  E/(l-y  ).  Substitution  of  y  and  dy  from  Equation  (50)  into 
Equation  (46),  ond  integration  over  the  holf-thickness  of  the  sheet,  gives  the 


rasult 


- /»4  <^4  , 

v2*t+/ 

)  J 

(51) 


^4 


wh«r«  Mp^  it  th«  monwnt  at  which  the  maximum  bending  stren  it  given  by 


(52) 
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and  0^  Is  the  maximum  bending  stress  corresponding  to  the  moment  M. 
Equations  (47)  and  (49)  evaluated  at  y  =  h/2  give  the  relation 


1 


(53) 


The  required  moment-curvature  relations,  from  which  the  effective  flexural  rigidity 
can  be  determined,  are  expressed  by  Equotions  (51)  to  (53). 

It  is  convenient  to  introduce  again  the  parameter  ^  , 
defined  by  Equation  (41),  With  this  relation  and  that  of  Equations  (23  )^ 

Equations  (51)  to  (53)  give  for  the  effective  flexural  rigidity. 


L 

I 


By  substituting  values  of  0\^/  •^PL  into  Equations  (54)  and  (55)  a  curve  of 
D'/D  vs  XM/Mq  can  be  constructued  for  use  with  the  incremental  procedure. 

This  result  may  be  compared  with  that  of  Equation  (43)  for  the  ideal  nonhordening 
mo^erial.  It  should  be  noted  that  the  ideal  nonhardening  material  con  be  con¬ 
sidered  as  0  limiting  case  of  Equation  (48)  as  n  approaches  infinity. 

For  the  onalysis  of  two-dimensional  convoluted  structures 
in  nondimensional  form  with  the  materiol  behovior  of  Equation  (48)  only  four  of 
the  parameters  introduced  need  be  specified.  These  are,  for  example,  O^p^/Ey  , 
Cp|_,  n,  and  A  ,  For  most  cases  of  interest  the  value  20  x  10*^  for  Cp|^  will 
probobly  suffice,  so  that  only  three  parameters  may  be  required.  Moreover,  for 
the  case  of  Equation  (43)  only  one  parameter  is  required. 

Examples  of  stress-strain  curves,  based  on  the  material 
behavior  of  Equation  (48),  are  shown  in  Figure  6  for  cases  of  strain-hardening 
exponents  corresponding  to  n  =  5,  n  =  10,  and  the  nonhardening  case,  n  =  ^  . 

The  other  material  porameters  assumed  ore  €p|^  -  20  x  10  ^  and  ~  10*^. 

The  corresponding  flexural  rigidity-moment  curves  are  shown  in  Figure  7. 

E.  Computational  Procedure 

The  procedure  by  which  the  incrementol  technique  is  used  to  compute 
the  deformation  behavior  of  a  convoluted  strip  is  described  in  the  following  steps: 

1.  The  typical  cantilever  section  is  divided  into  suitable 
increments  of  len^and  the  initial  values  of  Zj  and  0*  are  determined,  as  shown 
in  Figure  1c. 
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2.  A  small  increment  of  load  is  applied  and  the  corresponding 

deformations  ore  computed  from  small -deflection  elastic  theory.  For  simple  shapes 
classical  analytical  methods  can  be  used.  However,  if  the  shapes  are  more 
complicated  the  present  increment  technique  can  be  used.  For  this  purpose  the 
values  of  are  computed  from  Equation  (18),  taking  T  =  0  and  using 

the  initial  values  of  z.  for  z.  .  The  incremental  cantilever  slopes  and  deflections 
are  then  computed  from  Equations  (19),  (24),  and  (25),  and  the  resulting  deflections 

and  are  determined  from  Equations  (12)  to  (14). 

3.  The  deformations  corresponding  to  the  first  load  increment 

are  recomputed  using  the  results  of  Step  (2)  to  determine  better  values  for  the 
effective  incremental  moments,  The  values  of  T,  zj,  and  (dz/dT). 

are  computed  using  one-half  the  values  of  AT  and  •  If  the  nxiterial  is  in 

the  plastic  range  values  of  the  effective  flexural  rigidity  D!  are  computed  from  the 
average  moment  Tz*  and  the  appropriate  relations  of  Section  IV-D.  The  values  of 
’JTJfT  are  computed  from  Equations  (18)  or  (9),  depending  on  whether  the  material 
is  in  the  elastic  or  plastic  range.  The  deformations  are  then  computed  using 
Equations  (24)  and  (25)  or  Equotions  (32)  and  (33)  for  the  incremental  cantilever 
relations. 

4.  The  values  of  (dz/dT)*  are  computed  from  the  results  of 
Step  (3)  and  compared  with  the  initial  values  used  in  Step  (3).  If  the  agreement 
is  unsatisfoctory.  Step  (3)  is  repeated  using  the  final  values  of  (di/dT);  until 
ogreement  between  two  successive  values  is  attained. 
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5.  Subsequent  load  increments  are  applied  and  the  procedure  of 
Step  (4)  is  used  to  obtain  the  incremental  deformations. 

Although  an  iterative  technique  is  used  for  each  load  increment, 
convergence  appears  to  be  sufficiently  rapid  in  most  cases  that  the  method  con  be 
used  with  a  desk  calculator.  In  the  case  of  a  plastic  analysis  with  the  nonhardening 
material  behavior  some  difficulty  wos  experienced  due  to  the  severe  variation  of 
effective  flexurol  rigidity  with  moment,  as  illustrated  in  Figure  7.  However,  if 
the  load  increment  is  sufficiently  small  it  if  expected  that  this  difficulty  can  be 
avoided.  Moreover,  with  the  present  incremental  technique  it  is  just  as  easy  to 
use  Q  more  reolistic  moterial  behavior,  such  as  illustrated  in  the  other  curves  of 
Figures  6  and  7,  which  converge  more  readily.  A  more  complete  discussion  of  the 
computations  is  given  in  the  following  section. 

If  the  material  behavior  is  limited  to  the  elastic  range,  it  is  pasible 
to  use  a  direct  procedure  for  each  load  increment  and  to  correct  the  value  of  the 
load  increment  on  the  basis  of  energy  considerations.  For  this  purpose  it  is  con¬ 
venient  to  introduce  a  nxiteriol  efficiency  factor  defined  by 


elastic  strain  energy  in  structure 
elastic  Strain  energy  capacity 
corresponding  to  maximum  stress 


U 


where  U  is  the  elastic  strain  energy  per  unit  width  of  the  cantilever  beam  of 
Figure  1,  ^  is  the  total  length  of  the  beam,  and  ^  is  the  maximum  bending 
stress.  With  the  expression  for  U  in  terms  of  bending  moment,  given  by 


(56) 
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l;  ,  f  ^  =  -i-  r  Vr 


(57) 


ond  the  maximum  stress,  given  by 


g-  -  ^ 


where  is  the  maximum  bending  moment,  may  be  written 


(58) 


c/S  . 


(59) 


Since  the  elastic  strain  energy  in  the  structure  is  equal  to  the  work  done  during 
deformation. 


U=fT(r)J$'  , 


where  ^  is  the  total  horizontal  deflection,  as  indicated  in  Figure  Id. 
Equations  (56),  (58),  and  (60)  yield  the  relation 

/u^:  o 


3  > 

V  = . 

/Of  J 


(60) 


(61) 


The  maximum  bending  moment  is 

a  -  , 


(62) 
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where  is  the  maximum  moment  arm,  as  indicated  in  Figure  1c.  With  this 
result  Equations  (56)  and  (59)  can  be  written 


ond 


it 


(63) 


(64) 


It  will  be  shown  in  the  numerical  examples  that  V*  is  a  slowly 
varying  function  of  the  load  and  deflection.  Therefore,  by  determining  the  value 
of  a  given  load  increment  so  that  the  expressions  of  Equations  (63)  and  (64)  are 
equal,  a  more  accurate  result  con  be  obtained  on  the  basis  of  total  energy  absorbed. 
In  addition,  it  is  useful  to  evaluate  V'  as  a  function  of  load  and  deflection,  since 
this  quontity  is  a  measure  of  the  effectiveness  with  which  the  available  material  is 
being  utilized.  The  procedure  by  which  the  corrected  load  increment  is  determined 
is  described  below. 

After  application  of  a  particular  load  increment  1/^  is  determined 
From  Equotion  (63)  by  numerical  integration.  This  value  of  V  is  then  used  in 
Equotion  (64)  to  determine  the  corrected  lood  increment.  Let  4T^  dertote  the 
n*th  corrected  load  increment,  arMi  let  the  subscript  n  denote  values  of  the 
parameters  of  Equations  (60),  (63),  and  (64)  after  application  of  the  n-th  load 
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incramnt.  Then,  ossuming  a  linear  load-deflection  relation  for  each  load 
increment.  Equation  (64)  can  be  expressed  in  the  form 


(65) 


Th«  solution  for  A  from  Equation  (65)  is 


47; 


J 


where 


8h  s  a-7;.,  - 


X 


For  the  first  load  increment  corresponding  to  n  -  1,  the  solution  is 


477  = 


E±Jl  iii 
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Equotiom  (66)  and  (67)  may  b«  written  in  the  nondimensionol  form. 


(66) 


(67) 


(68) 


(69) 
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where 


F,  ComputatioT>al  Results 

The  accuracy  of  the  elemental  beam  approach  was  investigated  by 
the  application  of  this  part  of  the  method  to  the  small-deflection  analysis  of  Q 
semi-circular-arc  convolute  for  which  the  true  solution  is  known.  The  quarter 
circle  corresponding  to  the  cantilever  beam  of  Figure  1  was  divided  into  three, 
four,  and  five  increments,  and  the  resulting  load-deflection  relations  ore  compared 
with  the  known  solutions.  The  results  are  summarized  in  Table  I,  in  which  ^  arKl 
^  represent  the  horizontal  and  vertical  deflections,  respectively,  ar>d  A  is  the 
initial  rodius  of  curvature.  The  true  solutions  are. 


^4 
tttt  = 


W  " 


0.500 


(71) 


Table  I.  Small -Deflection  Solution  for  3,  4,  end  5  Beam  Increments 


No.  of 

%  error 

%  •rror 

Increments 

T7T5 

t/To 

in^ 

In  ^ 

3 

.757 

.505 

3.6 

1.0 

4 

.770 

.504 

1.9 

0.8 

5 

.775 

.501 

1.3 

0.2 
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It  is  seen  that  the  incremental  beam  technique  gives  rather  accurote 
results  for  relatively  few  increments.  For  example,  according  to  Table  I  the  use  of 
four  beam  elements  for  the  large-deflection  analysis  should  introduce  less  than  2% 
error  in  ^  and  less  than  1%  in  ^  ,  particularly  since  the  straight  beam  approxi¬ 
mation  improves  as  the  curvoture  decreases. 

Large-deflection  elastic  cases  were  carried  out  for  the  semi-circular- 
arc  configuration  and  two  U-shaped  configurations,  illustrated  in  Figure  8.  The 
principal  results  for  the  semi-circular-arc  are  shown  in  Figures  9  to  12,  and  those 
for  the  two  U-$hape$,  in  Figures  13  to  20.  The  results  include  the  load-deflection 
curve;  a  curve  of  ^/r  vs.  J^/r  which  indicates  the  "flattening”  of  the  convolute; 
a  curve  of  the  yield  parameter  (h/r)  /  (  ^^/Ey)  which  iruJicotes  the  deflection  at 
which  the  proportional  limit  is  reached;  and  the  variation  of  the  material  efficiency 
factor  V*  with  •^/r. 

For  comparison,  the  smol  I -deflection  solutions  for  the  U-shaped 
convolute  are  given  by 


r/n 


z  U/  ^  3  (a) 


•h  0.  ^7oe  -1 


(72) 


Plostic  COMS  w«r«  run  for  the  semi-circular-arc  configuration  with 
the  material  behavior  of  Figures  6  and  7,  and  two  values  of  the  yield  parameter  ^  . 
The  resulting  load-deflection  curves  are  shown  in  Figures  21  and  22.  As  mentioned 
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Figura  13.  Elotfie  iMwpad  Cotinguration, 
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Figur#  15.  Elpitic  Uflwyd  Gmfigurofion,  c^r  =  1 


Figm  15.  EloiHc  U"ihop<d  Configuration,  c^r  -  1 


Figure  17.  Elastic  U-shaped 


Figura  18.  Eloitic  U-shapad  Gmfiguration, 


1 


Figura  19  -  Elastic  U-shopod  OMifiguration, 
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Figwra  20  -  ElotHc  U-thopad  Connguratlon,  V*’  2 
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earlier,  some  difficulty  was  encountered  with  convergence  of  the  nonhardening 
material  case  due  to  the  sensitivity  of  effective  flexural  rigidity  with  bending 
moment.  The  curve  of  Figure  21  corresponding  to  n  =  oO  was  computed  with 
only  three  load  increments,  although  as  many  as  six  iterations  were  required  for  a 
particular  load  increment.  In  order  to  test  the  accuracy  of  these  calculations,  the 
same  case  was  repeated  using  load  increments  of  half  the  magnitude,  and  six  instead 
of  four  length  increments.  These  results  are  shown  in  Figure  23  and  indicate  only 
smell  deviations  from  the  original  solution.  The  curves  corresporKling  to  the  other 
cases  of  strain  hardening  were  computed  with  considerably  more  ease,  with  only 
one  or  two  iterations  being  adequate  for  some  of  the  load  increments. 

V.  THREE-DIMENSIONAL  CONVOLUTED  STRUCTURES 

In  a  two-dimen$ior>al  convoluted  structure,  such  as  the  strip  of  the  previous 
analyses,  deformation  arises  from  bending  in  the  plone  of  the  convolute  which  is 
essentially  constant  over  the  width  of  the  strip.  In  the  usual  three-dimensional 
convoluted  structure,  however,  additional  restraints  are  present  so  that  the  structure 
is  not  free  to  deform  by  o  uniform  bending  in  the  convoluted  plane.  For  example, 
in  the  common  bellows  of  Figure  24a,  extension  or  contraction  of  the  bellows  is 
occomponied  by  changes  in  the  diameter  which  produce  lateral  ''hoop"  stresses. 

These  stresses  impose  severe  restrictions  on  the  flexibility  of  the  bellows,  particularly 
in  the  cose  of  large  deflections.  Similarly,  in  the  deformotion  of  a  convoluted 
structure  from  a  cylindrical  to  o  conical  shape,  as  shown  in  Figure  24b,  the  non- 
uniform  circumferential  expansion  is  occomponied  by  twisting  of  the  longitudinal 
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This  produce  raitraints  which,  although  not  generolly  os  severe  os  In 
the  bellows,  con  significantly  affect  the  behavior  of  the  structure. 

For  the  present  analysis,  two  coses  of  the  circumferentially-convoluted 
cylinder  of  Figure  24b  are  treated.  The  first  case  represents  an  idealized  structure 
of  saw-tooth  configuration  which  expands  into  a  conical  shape,  and  illustrates  in 
a  simplified  monner  the  three-dimensional  restraining  effects.  The  second  cose 
deals  with  a  more  realistic  structure  in  which  the  convolute  has  one  of  the  U-shapes 
considered  earlier.  These  analyses  are  described  in  the  followir>g  sectior». 


A.  Circumferentiolly-Convoluted  Cylinder  with  Idealized  Sow-Tooth 
Configurotlon 

In  the  present  example  the  idealized  cylindrical  structure  is  restroined 
at  one  end  and  is  exparxJed  into  a  conical  shape  by  a  uniform  pressure  differential 
across  the  walls.  The  following  assumptions  ore  used  in  the  analysis; 

1 .  The  convolute  span  heigitis  small  compared  with  the  radius 

of  the  cylinder. 

2.  All  of  the  deformation  in  the  plane  of  the  convolute  occurs 
at  the  bends  so  that  only  the  angle  of  the  saw-tooth  is  changed. 

3.  A  linear  relation  exists  between  bending  moment  and 
convolute  (saw-tooth)  angle. 

4.  The  centerlirw  of  each  longitudinal  strip  remains  straight  so 
that  twisting  but  no  berxling  occurs  in  the  longitudinal  direction. 

5.  Small -deflection  torsion  theory  opplies  to  the  twisting  of  the 


strips. 


6.  The  material  behavior  is  elastic. 
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With  the  above  assumptions  expressions  can  be  derived  for  the  elastic 
energy  of  bending  and  tvristing  of  the  longitudinal  strips  and  the  work  done  on  the 
structure  during  expansion.  By  equating  the  work  done  with  the  total  elastic  energy 
the  relation  between  differential  pressure  and  cone  angle  con  be  determined. 

The  coordinates  and  some  of  the  geometric  parameters  are  shown  in 
Figure  25.  Because  of  symmetry  it  is  necessary  to  treat  only  one  of  the  longitudinal 
strips. 

The  twisting  energy  for  a  structure  comported  of  N  strips  is  given  by 

L 

0 


where  GJ  is  the  torsional  rigidity  of  the  strip.  From  Figure  25  we  have  the 
geometric  relations 

Z-n  R  —  N-l  ^ 


(74) 


and 


(75) 


Equotion  (75)  and  its  derivative  with  respect  to  x  give  the  result 


.  z 

i'/n  oC 


K,  -ZR,X, 


(76) 
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I  Figur*  25.  Coordinatm  and  Gaamatrle  Faromtlaft 

for  Convolutid  Cylindrieol  Strwetwrt 


Substitution  of  Equation  (76)  into  Equation  (73)  and  intogrotion  gives  the  result 


NdT  .  . 


(77) 


The  change  in  Uj  due  to  the  change  e(tk  in  the  cone  half-angle  is,  by 
differentiation  of  Equation  (77), 


olUy  ■m. 


Sixcf,  cp,) 


ZR, 


jL  O/M  <Pf  C/*t  ei 

R, 


(78) 


wh«r« 


■=  i/n"'  ^  sU  Cp,  J  . 


(79) 


The  bending  energy  for  the  strip  con  be  expressed  by 

L 

^  I  (80) 

c 


wh«r»  ^  it  Qn  •fftctivt  alottic  constont  for  the  joint.  With  the  substitution  of 
(p  from  Equation  (75)  into  Equation  (80),  the  result  of  the  integration  can  be 
expressed  by 


— —  f  u  (siit  'u)  ^  Zu  +■  2  iJf-u*-  Ci»'u 

Stii  ifh  sC  L 


Z<P,  (u  s'l^u  )J 

Si¥l^ 


(81) 
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Th*  ehangt  in  due  to  the  change  in  the  cone  half-angle  can  be  written  at 

a(  Sl^  ef^  of^ 


(82) 


Differentiation  of  Equation  (81)  with  Equations  (79)  and  (82)  gives  the  result 


(83) 


The  work  done  on  the  structure  by  the  differentiol  pressure  p  during 
the  expansion  is,  from  Figure  26, 

L 

^  ^  P  ^  (ZvRo/k)Xc^cC.  ,  (W) 

Integration  of  Equation  (84)  with  the  relation, 

R  ^  i-  X  SM  ^  ^  (85) 

gives  the  result 

^  VV  «  TT  /$,  l'p  (  /  f  ^  s/n  et  )  c/,t  .  (86) 

The  work  dorte  on  the  structure  can  be  equated  with  the  total  energy 
absorbed  in  the  structure,  to  give 

tf/lV  -  dUf^  ,  (87) 
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Substitution  of  Equations  (78),  (83)  and  (86)  into  Equation  (87)  gives  the 
required  preuure  relation, 

^^1  ^  f*  Sinef,  coiU 

NCT 


—  ^*^^**^  r. 


(88) 


This  example  illustrates  the  manner  in  which  three-dimensionol 
restraints  can  give  rise  to  twisting  as  well  os  bending  energy  absorption  by  the 
structure,  and  might  provide  o  first-approximation  solution  for  certain  convoluted 
configurations.  However,  o  more  realistic  onolysis  is  presented  in  the  followir>g 
section. 


B.  CircumferentiallyConvoluted  Cylinder  with  U-Shaped  Configuration 
The  present  arxilysis  is  similar  to  the  preceding  one  except  that 
lateral  berxiing  of  the  longitudinal  strips  is  treated  more  properly,  orKJ  the  resulting 
effects  on  the  twisting  of  the  strips  is  token  into  account.  A  typical  strip  of  the 
present  structure  has  an  S'-shaped  cross  section  which  deforms  as  shown  in  Figure  27. 
The  lateral  deformation  of  the  strip  at  a  particulor  axial  location  is  determined  from 
the  two-dimensional  analysis.  It  is  seen  from  the  sketch  of  Figure  27  that  the  angle 
of  twist  varies  over  the  cross  section  of  the  strip,  os  compared  with  the  uniform 
twisting  in  the  previous  analysis.  This  variable  twisting  is  taken  into  account  by 
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trvoHng  th«  strip  os  though  it  woro  mode  up  of  many  smaller  strips,  each  of  which 
twists  a  different  amount.  From  symmetry  it  is  necessary  to  treat  only  half  of  the 
S-shaped  strip,  os  in  the  two-dimensional  onalysis. 

The  following  assumptions  are  made: 

1 .  The  convolute  span  height  is  small  compared  with  the 
radius  of  the  cylinder. 

2.  The  center  line  of  each  longitudinal  S-shaped  strip  remains 
straight  os  the  cylinder  deforms  into  a  conical  shape. 

3.  The  load-deformation  behavior  due  to  lateral  bending  at 
any  axial  location  is  the  same  as  that  determined  in  the  two-dimensionol  analysis. 

4.  Smoll-deflection  torsion  theory  applies  to  the  twisting  of 
each  element  of  the  S-shaped  strip,  ar>d  a  constant  torsional  rigidity  relotes 
twisting  moment  per  unit  length  of  cross  section  with  the  local  angle  of  twist. 

5.  The  material  behavior  is  elastic. 

With  the  above  assumptions  the  twisting  energy  of  the  structure 
per  unit  length  in  the  axial  direction  can  be  written 


wh«r«  N  !i  the  number  of  S-thaped  strips  and  GJ  is  th«  torsional  rigidity  of  the 
strip  per  unit  length  of  cross  section.  The  angle  of  twist  ^(x,s)  is  a  function  of 
both  the  axial  coordinate  x  and  the  cross  section  coordinate  s.  For  o  strip  of 
narrow  crou  section  the  torsiortal  rigidity  per  unit  length  is,  from  Reference  D, 


(yj  - 


(90) 
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If  It  it  auuimd  tliat  th«  functional  form  of  ^  (x,i)  with  s  it  tho  tamo  for  oil 
valuot  of  X,  Equations  (89)  and  (90)  con  bo  written 


[  of ‘Pm  , 

lu - 3 -  1  > 


(91) 


whtra  ky  is  ind«p«nd«nt  of  x  orKl  is  doftnod  by 


<"«'  t  is  tho  maximum  ongle  of  twist,  corresponding  to  s  -  .  Although 

oppeors  to  be  a  function  of  x  by  Equation  (92),  it  will  be  shown  in  the 
example  that  it  is  nearly  invariant  in  x. 

Integration  of  Equation  (91)  from  0  to  L  gives  for  the  total 
twisting  energy. 


/ 

f  r  ^  7  ^  y  /  A  I 

3L  J 


expansion 


For  the  U-shaped  configuration  discussed  earlier,  the  fractional 
is  related  to  the  present  structural  parameters  by 

f/A.  =  ^  Ct'otoC  ^ 


(92) 


(93) 


(94) 


and 

=  k,  ‘^(t)  ,  <*» 
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so  that  Equation  (93)  can  ba  written. 


U, 


A/  (y  ^  ^  ,  (*  r<^e(l,  -7  i  „ 

=  "7^""^  15^)  J 


(96) 


Th#  functional  ralotion  for  d  ^/d(^/^)  vs.  ^ can  be  approximated  by 

^  c 


(97) 


where  a^  ond  b^  are  determined  from  the  two-dimensional  analysis.  Substitution 
of  Equation  (97)  into  Equation  (96)  and  differentiation  with  respect  to  gives 

S R,  L'  (/i,) +  3  h  (m,)  5'*»o(jco<oi . 


This  expreuion  can  be  used  with  the  differential  energy  relation,  as  in  the  preceding 
anolysis,  or  integrated  to  give  the  total  twisting  er^ergy  as  a  furKtion  of  ^  , 

/I  t  /  ^  7  I  z.  j  L  ^  c-'N. 


e>0. 


The  bending  energy  con  be  expressed  in  the  form. 


(100) 
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wh«r«  T  is  tha  circumftrtntial  toad  per  unit  length  and  is  the  reference  load 
of  Equation  (23).  The  function  T/T^  {%/K  )  can  be  approximated  from  the  two- 
dimensionol  analysis  by 


(101) 


Substitution  of  Equation  (101)  into  Equation  (100)  ond  integration  gives  the  result 


dK 


which,  with  Equation  (94),  can  bo  written 


cl  Ug 


(102) 


(103) 


Integration  of  thii  aquation  from  ^/A.  =  0  to  ^ /k  -  (l/R|)sino^  gives 
the  result 


-  3  A/4  ^,7;  y.  ^  ^  ^ 


(104) 


The  differential  bending  energy  expression,  analogous  to  Equotion  (98),  can  then 


be  written 


(105) 


With  th«  •xpraiflon  for  from  Equation  (23),  and  tho  rolotioni. 


g 


(106) 


Equationi  (99)  and  (100)  can  bt  combined  to  give  the  total  energy  input  W  as  a 
function  of  / 


W  =  (^r  + 


(107) 


+4n-y)  + 


Tho  proMnt  mothod  is  illustrotod  with  on  oxompio  bosod  on  tho 
U>shapod  configuration  with  q/^  *  1,  mod  in  tho  Iwo-dimonsiortal  analysis. 

Tho  dotormination  of  tho  constants  kj,  a^,  b^,  02,  ond  ^>2  from  tho  two-dimonsional 
analysis  is  doseribod  in  Appendix  A.  Those  results  are  os  follows: 


k^  =  0.460,  o,  =  0.647,  b,  =  0.121, 

a  -  0.213,  b.  »  0.0612 
2  2 


(108) 


With  those  values,  and  Jt/ji  *  2.571  for  this  configuration.  Equation  (107)  can 
be  written 
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If  th«  span  haight  of  tho  convolutod  structure  is  not  small  compared 
with  the  radius  of  the  cyltrMler,  os  ouumed  in  the  foregoing  anolysis,  the  resulting 
equotions  are  r^ot  strictly  applicable,  due  to  the  significant  curvature  imposed  on 
the  (gross)  structure.  This  curvoture  gives  rise  to  additional  twisting  deformations 
os  the  cylir^der  is  expanded  into  the  conical  shape  and  also  effects  the  berxJing  in 
the  plones  of  the  convolute.  However,  orte  simple  correction  con  be  applied  in 
order  to  improve  the  theory  for  this  case. 

Figure  28  illustrates  the  monner  in  vrhich  a  circumferentiol  strip  is 
deformed  durirtg  forming  of  the  cylirKler  from  a  straight  convoluted  strip.  It  can  be 
seen  thot  sif^ificont  bending  occurs  during  the  forming  operation.  Here,  "A** 
represents  the  ur>deformed  shape  and  "B”  is  the  shope  in  the  cylindrical  structure. 
If  the  cytindricol  structure  is  not  stress-relieved  prior  to  expansion,  as  in  the 
experiment  described  in  the  following  sections,  some  of  the  elastic  strain  energy 
will  be  releosed  during  expansion.  Conversely,  if  the  structure  is  stress-relieved. 
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lonw  additional  energy  will  be  required  to  expand  it. 
the  strip  of  Figure  28  from  "A"  to  “B"  it 


The  energy  required  to  deform 


(110) 


where  iZ  is  the  length  of  the  S-thaped  strip  and  Up  is  the  total  energy  required 
to  form  the  cylinder.  Since  the  total  angle  of  rotation  of  the  ends  of  the  strip  is 
2'ir  /N, 


(in) 


and  Equation  (110)  becomes 


NX 


(112) 


where  the  relation  of  Equation  (3)  has  also  been  used.  During  expansion  of  the 
cylirder  Into  the  cone  the  angle  of  rotation  of  the  ends  of  eoch  S-shaped  strip  is  reduced 
by  the  foctor  cot cL  ,  where  o6  is  the  cone  half-angle.  *  The  change  in  Up 
during  expansion  is,  by  Equations  (1 10)  and  (112),  proportionol  to  the  squore  of  the 
angle  of  rotation.  Thus,  the  releose  in  energy  for  the  nonrelieved  case  is  given  by 


AUff  Uf:  (l^  )  , 


(113) 


*  Dut  to  tha  twisting,  th*  strip  is  not  ratatad  uniformly  ovar  its  langth.  Hovravar, 
sine*  this  is  a  socortd  ordar  affect,  tha  avaraga  cona  ongla  is  tokan  for  simplicity. 
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or,  wUh  Equation  (112), 


AUp 


ir^ 

7T 


A/  X 


.  X 
S/*I 


(114) 


Subtraction  of  this  rasult  from  Equation  (107)  yields 


(115) 


For  the  U-shaped  configuration  with  a/A,  -  1  and  the  parameters  of  Equations  (108), 


Equation  (109)  becomes 


w  = 


0.005*92 


\*- 

a/2-  L  l) 


i-/3. 


0.0 /li 


(r,) 


(116) 
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If  there  is  an  <qi>preciable  difference  between  the  inner  and  outer  radii  of  the  cylinder, 
it  is  expected  that  Equations  (115)  and  (116)  will  still  be  only  approximately  correct 
due  to  the  twisting  effects  mentioned  earlier.  As  the  cylinder  expands  into  a  conical 
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thap«  the  nonuniform  circumferential  elongation  produces  a  nonuniform  flattening 
of  the  convolute,  with  the  maximum  flattening  occurring  at  the  expended  end. 
Consequently,  the  inside  generators  of  the  cylinder  rotate  through  o  greater  angle 
than  those  on  the  outside,  producing  additional  twisting  deformations.  Thus,  the 
results  of  the  two'-dimensional  bending  analysis  are  not  directly  applicable  since 
some  of  the  circumferential  elongation  is  due  to  this  twisting  deformation.  A  more 
proper  analysis  of  the  problem  should  treat  the  bending  and  twisting  deformations 
simultaneously. 

VI.  EXPERIMENTS 

Tests  were  conducted  with  several  two-dimensiorKil  specimens  and  a  three- 
dimensional  cylindrical  structure  similar  to  that  described  in  the  preceding  section. 
Descriptions  of  the  test  procedures  or>d  the  experimentol  results  are  given  in  the 
following  sectiortt. 

A.  Two-Dimensional  Strips 

Lood-deflection  tests  were  cor>ducted  on  two-dimensional  Type  304 
stainless  steel  specimens  with  semi-circulor-arc  and  U-shaped  configurations,  os 
used  in  the  previous  numerical  examples.  Samples  of  the  specimens  ore  shown  in 
Figure  29.  The  specimens  were  fabricoted  from  0.002"  and  0.003"  onneoled 
sheet  materiol.  Some  of  the  specimera  were  tested  in  the  as-formed  state  arKi  others 
were  procen  onneoled  prior  to  testing.  The  rodius-to-thickness  rotios  were  selected 
to  give  approximately  20  percent  expansion  prior  to  yielding  for  the  semi-circulor-orc 
specimens  otkI  63  percent  exponsion  for  the  U-shoped  specimens.  A  summary  of  the 
specimens  tested  is  given  in  Table  II. 
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Figurt  29.  Two-Dimensional  Specimens  and  Small  Weights 


Table  II.  Two-Dimensional  Specimens 


Specimen  No. 

Type 

Nominal  Thickness 
(in.) 

Condition 

NomirKil 

length/radius 

1 

U 

.003 

as  formed 

4 

2 

U 

.003 

os  formed 

4 

3 

U 

.003 

annealed 

4 

5 

semi-circle 

.003 

as  formed 

8 

6 

semi-circle 

.003 

as  formed 

8 

8 

semi-circle 

.003 

annealed 

8 

9 

U 

.002 

as  formed 

4 

10 

U 

.002 

as  formed 

4 

12 

U 

.002 

onnealed 

4 

13 

semi-circle 

.002 

os  formed 

8 

14 

semi-circle 

.002 

as  formed 

8 

15 

semi-circle 

.002 

anrieoled 

8 

Width  =  4.00  inches  (all  specimens) 

Nominal  (thickness/radius)  -  0.00730  (all  specimens) 
EstiitKited  variation  in  thickness  -  +  0.00005  (all  specimens) 
Estimated  variation  in  radius  =  ^  10%  (as-formed  specimens) 

Annealed  specimens  exhibit  considerable  distortion. 
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The  specimeni  were  suspended  on  knife-edges,  os  shown  in  Figures  30 
and  31,  in  order  to  produce  zero  moment  loading  conditions  ot  the  ends  of  the  active 
length  of  the  specimens.  This  was  occomplished  through  fhe  use  of  fine  wires  wrapped 
around  thin  metal  strips  and  passed  through  small  holes  in  the  ends  of  the  specimens. 
Loading  was  accomplished  by  the  addition  of  calibrated  weights  to  a  light  pan 
suspended  from  the  lower  end  of  the  specimen.  Deflections  were  measured  with  a 
vernier  height  gage  and  attaching  telescope  by  sighting  on  the  upper  and  lower  ends 
of  the  specimen. 

B.  Three-Dimensional  Cylindrical  Structure 

A  circumferential  I  y-convoluted  cylindrical  structure  was  formed  by 
joining  the  ends  of  two  convoluted  strips.  The  material  and  dimensions  were  the  same 
as  U-shaped  specimens  Nos.  1  and  2,  and  the  total  length-to-radius  of  the  strips 
before  joining  was  40.  Manner  of  loading  is  illustrated  in  Figures  32  to  35.  The 
energy  required  to  deform  the  structure  into  a  conical  shape  was  measured  for 
various  cone  angles  ond  compared  with  the  theoretical  predictions. 

Loading  was  accomplished  by  the  use  of  ten  mass  balanced  frames 
secured  to  the  outermost  generators  of  the  cylindrical  structure  and  pivoted  about 
the  bose  in  rodial  plones.  Fine  wires,  used  to  attach  the  frames  to  the  structure, 
were  passed  over  slots  in  the  frames  and  joined  to  a  common  point  below,  from  which 
a  light  balance  pan  was  suspended*  Calibroted  weights  were  added  to  the  pan  and 
the  corresponding  inner  and  outer  diameters  of  the  upper  end  of  the  structure  were 
measured.  The  changes  in  diameters  were  used  to  compute  the  cone  angle  and  the 
movement  of  the  pan,  from  which  the  corresponding  energy  inputs  were  computed. 
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Figure  33.  Convoluted  Cylindrical  Structure  Loaded  to  Partial  Expaiuion 
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34.  Convoluted  Cylindrical  Structure  View  of  Loading  Segments 
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Figure  35.  Convoluted  Cylindrical  Structural  Test  Setup,  Side  View 
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C.  Expertmantql  Rtsulti 


1,  TwO">Difnantionql  Spaciment 

Tha  load-daflaction  curves  for  tha  two-dimansional  spacimens 
ora  shown  in  Figuras  36  to  39.  Tha  corresponding  theoretical  curves  ore  also  shown, 
bosad  on  tha  nominol  dimensions  and  values  of  29.4  x  10^  psi  for  E  and  0.29  for 
V  •  Tha  considerable  scatter  in  tha  data  can  be  attributed,  in  part,  to  variations 
in  the  dimensions  of  the  specimens.  This  is  suggested  by  the  fact  that  the  initial 
slopes  of  the  curves  ore  also  in  disagreement  with  the  small -deflection  values  and, 
consequently,  these  curves  appear  to  be  rotated  about  the  origin.  In  order  to 
obtain  more  valid  comparisons  with  the  large-deflection  theory,  the  data  hove  been 
replotted  in  nondimensional  form  by  using  the  initial  slopes,  with  the  small -deflection 
theory  solutions,  to  compute  values  of  the  reference  load  T^.  These  curves  ore 
shown  in  Figures  40  to  43. 

2.  Convoluted  Cylindrical  Structure 

Figure  44  shows  the  energy-deformation  curve  for  the 
convoluted  cylindrical  structure,  in  which  energy  input  is  plotted  against  cone  half¬ 
angle.  The  corresponding  theoretical  curve,  based  on  Equation  (115),  is  also  shown. 
For  these  computations  the  nominal  dimensions  of  the  two-dimensional  strip  were  used, 
with  a  measured  volue  of  3.142  inches  for  the  average  initial  cylinder  radius,  R^« 

VII.  CORRELATION  OF  THEORY  AND  EXPERIMENT 
A.  Two-Dimensional  Analysis 

The  scatter  in  the  curves  of  Figures  36  to  39  con  be  attributed,  in  part, 
to  the  variation  in  specimen  dimensions,  as  noted  above.  For  example,  the  10% 


-80- 


ri&.  37  LOAP- DEP lection  CURVES  FOR  0.002- INCH  O  -  5HA.PED  SPECIN^EN'S 


Fl6,39  UOAD- DEFLECTION  CuRvES  FOR  0.002-INCH  SEM  I- C  IRCULAR  -  ARC  SPECINACNS 


hH 


ri6,  41  NONDIMCnSIONAU  load-oeflection  curves 

FOR  0.003  -  INCH  CIRCULAR  ARC  SPECIMENS 


FIG.  42  NONDIME NSIONAL  LOAD-DEFLECTION  CURVES 
FOR  0.002-INCH  U-SHAPED  SPECIMENS 


r 
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FIG.  4  3  NONOIMENSIONAL  LOAO-DEFLECTION  CURVES 
FOR  0.002-INCH  SEMI-CIRCULAR  ARC  SPECIMENS 
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variation  in  radius  can  produce  over  20%  variation  in  load,  and  the  0.05  mil 
variation  in  thickness  can  give  rise  to  variations  of  approximately  5%  and  6% 
in  load,  respectively,  for  the  2  and  3  mil  specimens.  These  variations  can 
accumulate  and  thereby  account  for  most  of  the  scatter  in  the  date. 

A  considerable  improvement  results  in  the  "rotated"  curves  of 
Figures  40  to  43.  It  is  interesting  to  note  the  similarity  in  the  curves  for  the  3  mil 
coses  and  the  closer  correlation  than  in  the  2  mil  cases.  This  would  be  expected 
since  the  dimensional  tolerances  that  can  be  attained  are  about  the  same  for  both 
cases,  but  a  given  dimensional  error  in  a  2  mil  case  has  a  greater  effect  than  in  a 
3  mil  case. 

It  is  expected  that  the  consistently  lower  curves  for  the  as*formed 
specimens,  compared  with  the  anneoled  specimens,  are  due  to  the  Bouschinger 
effect.  *  This  effect  is  related  to  peculiarities  in  the  stress-strain  curve,  particularly 
the  earlier  yielding  and  lowering  of  the  curve,  when  a  polycrystalline  metal  is 
plastically  deformed  in  one  direction  ond  then  looded  in  the  reverse  direction.  This 
situation  occurs  with  the  convoluted  specimens.  The  forming  operation  produces 
plostic  bending  deformation  in  one  direction,  and  elongation  of  the  specimens 
results  in  bending  in  the  reverse  direction.  Moreover,  a  mild  annealing  is  known 
to  remove  these  effects  and  restore  the  metal  to  its  normal  behavior. 

The  disogreement  between  the  theory  and  experiment  that  persists  in 
the  "rototed"  curves  of  Figures  40  to  43,  besides  possible  Bauschinger  effects  in  the 


*  See,  for  exomple.  Reference  B,  Chapter  12,  and  Reference  E,  Section  3.4* 
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Qs-formed  cases,  can  probably  be  attributed  to  dimensional  nonuniformities. 

E.g.,  the  curves  were  not  truly  circular,  and  variations  in  curvatures  and  other 
dimensions  occurred  within  a  single  specimen.  This  is  particularly  true  of  the 
annealed  specimens  in  which  considerable  distortion  resulted  from  the  annealing 
process.  Nevertheless,  the  overall  ogreement  appears  to  be  quite  good  so  that 
the  theory  appears  to  be  substantiated,  at  least  in  the  elastic  range. 

It  should  be  mentioned  that  some  plastic  yielding  occurred  over  the 
latter  portion  of  the  lood-'deflection  curves.  No  attempt  was  made  to  compute  this 
effect  with  the  exception  of  the  plastic  computations  on  the  semi-circular-arc 
configuration  presented  in  Section  IV-F.  The  lack  of  adequate  stress-strain  data 
for  very  thin  sheet  material  and  the  scatter  in  the  data,  coupled  with  the  slight 
effect  yielding  has  in  this  range,  would  make  correlations  quite  difficult.  The 
effects  of  plastic  deformation  can  be  evaluated  more  easily  with  larger  specimens 
in  which  dimensions  and  material  properties  can  be  more  closely  controlled. 

Convoluted  Cylindrical  Structure 

The  curves  of  Figure  44  disogree  considerably  over  the  lower  range 
of  expansion  but  show  rather  good  agreement  at  the  higher  range.  The  maximum 
error  is  about  65%  and  the  average  error  is  about  30%.  This  error  can  be  attributed 
to  the  appreciable  difference  between  the  inner  and  outer  radii  of  the  cylinder  and 
the  consequent  inadequocy  of  the  theory  in  this  case.  As  mentioned  in  Section  V-B, 
this  difference  in  radii  gives  rise  to  additional  twisting  effects  which  are  not  accounted 
for  in  the  present  analysis.  A  more  proper  armlysis  of  the  deformotion  should  treat  the 
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large-deflection  bending  and  twisting  of  the  circumferential ly-convoluted  strips 
simultaneously. 

A  more  quantitative  substantiation  of  the  source  of  error  arises  from 
a  consideration  of  Equation  (115)  with  the  dimensions  of  the  structure  and  the  basic 
assumption  that  the  difference  in  radii  is  small  compared  with  the  overage  radius. 
According  to  Equation  (115)  the  energy  is  approximately  inversely  proportional  to 
the  square  of  the  average  radius,  .  In  the  actual  structure  the  inside,  outside, 
and  average  radii  were  2,34  inches,  3.94  inches,  and  3.14  inches,  respectively. 
Thus,  if  the  inner  and  outer  radii  were  used  to  compute  the  energy,  rather  then  the 
average  radius,  variations  of  approximately  80%  and  -37%,  respectively,  would 
result.  In  view  of  this  it  appears  that  the  theory  gives  a  fairly  good  first  opproxi- 
motion,  particularly  for  the  upper  range  of  expansions.  Moreover,  it  would  appear 
that  considerable  improvement  in  the  theory  would  result  os  the  difference  in  radii 
were  decreased* 

VIII .  CONCLUSIONS  AND  RECOMMENDATIONS 

Theoretical  methods  of  anolysis  hove  been  developed  for  two-dimensionol 
convoluted-type  structures  subjected  to  large  deformations,  with  material  behavior 
in  both  the  elastic  and  plastic  ranges.  Computations  have  been  carried  out  for 
structures  of  several  typical  configurotions,  utilizing  numerical  techniques  suitoble 
for  use  with  a  desk  colculator.  Results  of  the  two-dimensional  analysis  have  been 
applied  to  the  analysis  of  a  typical  three-dimensional  convoluted  cylindricol 
structure. 
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Ttsts  have  been  performed  on  two-dimensional  specimens  and  a  three- 
dimensional  cylindrical  structure  similar  to  thot  used  in  the  analysis.  The  testeappear 
to  substantiate  the  two-dimensional  analysis  within  the  elastic  range  of  behavior. 

The  results  with  the  convoluted  cylindrical  structure  indicate  that  the  simplified 
theory  developed  can  be  used  for  a  first-approximation  analysis  of  a  considerably 
more  complex  structure.  Thcyfurther  indicate  that  a  considerably  more  refined 
arKilysis  will  be  required  in  order  to  accurately  predict  the  behavior  of  a  three- 
dimensior>al  structure  of  the  type  tested. 

In  view  of  the  limited  experimentation  that  could  be  performed  under  the 
present  study  the  following  investigations  are  recommended  in  order  to  amplify  and 
extend  the  present  findings: 

A.  Tests  with  larger  two-dimensionol  specimens  so  that  dimensions  and 
material  behavior  can  be  more  carefully  controlled 

B.  Tensile  tests  to  determine  stress-strain  behavior  of  specimen  sheet 
moterials 

C.  Simple  bend  tests  for  independent  determinations  of  plastic  bending 
behavior 

D.  Compression  tests  with  as-formed  specimens  to  determine  initial 
slopes  and  evaluate  Bouschinger  effects 

E.  Tests  with  the  some  specimen  in  the  as-formed  and  annealed 
conditions  in  order  to  evaluote  Basuchinger  effects 

F.  Refined  analyses  of  the  convoluted  cylindrical  structure  tested  in 
order  to  account  for  the  large  difference  in  radii 

G.  Tests  with  convoluted  cylirtdrical  structures  having  smaller  differences 
in  the  radii  in  order  to  determine  the  range  of  applicability  of  the  simplified  theory 
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In  addition  to  the  above  investigations  it  is  recommended  that  future 
theoretical  and  experimental  studies  be  carried  out  on  three-dimensional  convoluted 
cylindrical  and  conical  structures  in  which  (a)  bending  of  the  structure  occurs  in 
the  axial  as  well  as  in  the  circumferential  directions;  and  (b)  the  convolute  span 
height  varies  axially  so  as  to  produce  preferential  expansions  of  the  structure.  In 
particular,  attempts  should  be  made  to  determine  if  a  simplified  theory  based  on  the 
two-dimensional  structure  can  be  developed,  os  in  the  present  study. 

It  should  be  noted  that  the  development  of  a  successful  theory  for  analyzing 
the  structures  described  by  (a)  ond  (b)  would  constitute  a  significant  step  toward 
understanding  the  general  behavior  of  convoluted  shell  structures.  As  seen  from 
the  simplified  analysis  of  the  convoluted  cylindrical  structure  in  the  pre.^ent  study, 
the  general  treatment  of  convoluted  shell  structures  will  be  extremely  difficult. 
Therefore,  if  it  can  be  developed,  a  simplified  theory  based  on  an  extension  of  the 
present  study  would  be  quite  useful  for  design  purposes. 
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APPENDIX  A 


DttTminotion  of  Pqrometers  for  Convoluted  Cylindricql  Structure 
from  Two*"bifnensional  An^yTis 


The  determination  of  the  parameters  a^,  02/  and  b2/  defined  in 

Equations  (92),  (97),  and  (101),  is  described  below.  The  porameter  k^  is  a 
measure  of  the  distribution  of  rotation  over  the  cross  section  of  a  typical  longitudinal 
strip  of  the  cylinder  and  is  determined  from  the  changes  in  slopes  of  the  two-dimensional 
structure,  as  illustrated  in  Figures  1  and  27.  The  parameters  and  b^  relate  the 
maximum  rotation  of  the  strip  to  the  fractional  elongation,  and  a2  ond  b2  describe 
the  two-dimensional  load-deflection  curve. 

It  wos  assumed  in  Section  V-Bthat  k^,  as  defined  in  Equation  (92)  is,  in 
effect,  independent  of  the  circumferential  elongation  of  the  structure.  The  validity 
of  this  ouumption  for  the  U-shaped  configuration  used  in  the  numerical  example  and 
in  the  experiment  is  demonstrated  by  Figures  45  and  46.  Figure  45  shows  distributions 
of  the  changes  in  slope,  ^  ,  over  the  cross-section  of  the  typical  section  for  several 
elongations,  as  determined  from  the  two-dimensionol  arwilysis.  Figure  46  shows  the 
correspor>ding  distributions  of  the  incremental  changes  in  slopes,  normalized  on  the 
bosis  of  the  maximum  values.  The  slight  variation  in  the  curves  of  Figure  46  over  the 
range  of  expansion  shown  justifies  the  inltiol  assumption.  The  value  for  k^  was 
obtained  from  the  data  of  Figure  46  and  Equation  (92),  which  con  be  approximated 
as  follows: 


(A-1) 
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The  average  values  of  i<P/^  cP  were  used  in  Equation  (A-1)  and  was 

m  I 

determined  by  numerical  integration. 

The  parameters  a^  b^,  and  q^,  b^  were  obtained  by  fitting  the  curves 

of  <P  vs.  and  T/T  vs.  J^//t ,  determined  from  the  two-dimensional 

m  o 

analysis,  to  Equotions  (97)  and  (101),  respectively.  The  original  curves  and  fitted 
expressions  are  shown  in  Figures  47  and  48. 
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